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I.  Project  Objective 


The  objective  of  this  program  was  to  develop  a  new  approach  to  modeling  terrain 
using  mathematical  concepts  from  fractal  and  attractor  theories.  Key  parameters  of  the 
model  were  extracted  from  fingerprint  analysis  based  on  scale  space  smoothings  of 
elevation  profiles. 


II.  Scale  Space  and  Fingerprint  Background 

Description  of  terrain  involves  identification,  over  a  broad  range  of  physical  scales,  of 
the  spatial  extent,  prevalence,  and  relative  importance  of  elements  such  as  peaks,  ridges, 
and  valleys.  These  elements  are  delimited  by  fundamental  edges  of  the  topography,  and 
so  it  is  natural  to  propose  terrain  modeling  based  on  a  scale  space  filtering  of  elevation 
data.  A  scale  space  image  is  a  set  of  progressively  smoothed  versions  of  the  original 
topographic  array:  as  the  smoothing  scale  increases,  features  of  the  terrain  disappear  until 
only  a  dominant  terrain  shape  remains.  A  plot  of  the  points  of  inflection,  or  edges,  within 
the  scale  space  image  results  in  a  lingerprint’  of  the  terrain.  The  fingerprint  locates  points 
of  inflection  of  the  original  data,  provider-  a  i/?ntext  in  which  the  inflection  points  are  paired 
and  assigned  to  a  topographic  feature,  relates  topographic  features  of  the  curve  to 
each  other. 

Witkin  developed  the  scale  space  fingerprint  techn’que  in  order  to  relate  signal 
descriptions  across  different  filtering  scales  (Ref.  2).  Its  mathematical  foundation  lies  in  the 
theory  of  heat  flow.  A  signal  f  (in  our  case  f  is  a  slice  of  a  topographic  elevation  array)  can 
be  regarded  as  an  initial  temperature  distribution  on  the  real  line.  As  heat  flows,  the 
temperature  distribution  smoothes;  weak  features  (local  hot  or  cold  spots)  vanish  in  a 
short  time  while  stronger  features  persist  for  longer  times.  A  feature  is  delimited  by  its  left 
and  right  inflection  points;  when  these  points  merge,  the  feature  vanishes.  No  new 
features  (sources  or  sinks)  arise  during  heat  flow.  The  complete  temperature  record  is 


1 


called  the  scale  space  image,  and  the  graph  of  inflection  points  over  time  is  called  the 
fingerprint.  Scale  space  techniques  have  recently  been  applied  to  edge  detection  and 
pattern  recognition,  to  the  cartographic  matching  of  shorelines,  to  histogram  analysis  of 
normal  mixtures,  and  to  modeling  and  analysis  of  hyperspectral  data  curves. 

In  applying  scale  space  techniques  to  terrain  analysis,  the  input  signal  f(x)  is  an 
elevation  profile  obtained  by  either  an  east-west  or  a  north-south  scan  from  7.5  minute 
United  States  Geological  Survey  Digital  Elevation  Model  (DEM)  records.  One  x-unit 
corresponds  to  30  meters.  G°  is  a  mean  zero  Gaussian  mask  with  standard  deviation 
a  >  0: 


G%)  =  J_  e 
o 

When  we  convolve  f  with  G°,  we  obtain  the  scale  space  image 

f  oo 

F(x,o)  =  (f*G°)(x)  =  f{x-y)  GP(y)  dy 

For  each  fixed  a,  F(x,a)  is  a  smoothing  of  f(x).  An  edge  occurs  where  a  function  is 
changing  most  rapidly.  Equivalently,  edges  are  points  of  inflection,  where  the  slope  has  a 
local  maximum  or  minimum.  The  fingerprint  tracks  the  location  of  inflection  points  of  F(x,a) 
against  smoothing  scale  a;  thus,  it  is  the  plot  of  all  points  (x,cy)  where  changes 

sign  (Fig.  1). 

The  input  f  of  our  terrain  analysis  consists  of  discrete  elevation  data,  yet  our  scale 
space  analysis  involves  differentiation  to  get  slopes  and  inflection  points.  Gaussian 
convolution  efficiently  deals  with  this  problem,  since  the  result  of  convolution  of  a 
bounded  measurable  f  with  a  smooth  kernel  such  as  the  Gaussian  produces  an  output 
function  that  inherits  the  differential  properties  of  the  kernel.  In  particular, 

^  <(y)G''(x-y)dy  =  J.~  ♦(y) 
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Fig.  1.  An  elevation  profile  (top)  and  corresponding  fingerprint  (bottom)  for  a  9,000  m 
trace  from  the  Lassen  Volcanic  Wilderness  region  of  California.  Location  of  features  with 
closing  scale  o^  7  are  superimposed  on  the  elevation  profile.  Each  unit  of  a  corresponds 
to  a  data  spacing  unit  of  30  meters. 
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and  if  L  is  any  linear  differential  operator,  then 


L(f*G®)  =  f*LG° 

Since  we  use  numerical  convolution  based  at  each  data  point,  we  are  able  to  calculate 
derivatives  at  each  data  point.  Gaussian  convolution  is  particularly  well  suited  to  handling 
differentiation  of  discrete  imagery,  since  in  a  single  operation,  such  as  f  *  (G°)',  we  both 
smooth  and  differentiate.  The  degree  of  smoothing,  or  filtering,  is  controlled  by  the 
standard  deviation  parameter  a  of  the  Gaussian  kernel. 

The  Gaussian  G®  shares  a  variation  diminishing  property  with  other  diffusion  kernels; 
the  number  of  zero-crossings  of  f  *G°  is  not  more  than  the  number  of  zero-crossings  of  f. 
Consequently,  the  number  of  inflection  points  of  F(x,ct)  also  cannot  increase  when  a 
increases.  This  last  property  imparts  a  characteristic  appearance  to  a  fingerprint. 
Qualitatively,  a  fingerprint  appears  as  an  assemblage  of  closed  arches  (Fig.1,  Fig.  2).  Each 
arch  is  composed  of  a  pair  of  zero-crossing  contours  of  one  solid  (darker) 

contour  and  one  dotted  (lighter)  contour.  As  we  cross  a  solid  contour  from  left  to  right,  the 
smoothed  curve  F(x  ,o)  changes  from  concave  down  to  concave  up.  As  we  cross  a  dotted 
contour  from  left  to  right,  F(x,a)  changes  from  concave  up  to  concave  down.  As  the 
smoothing  scale  a  increases,  zero  crossing  contours  pair  off,  merge  to  form  an  arch,  and 
die  out  at  the  apex  of  the  arch.  The  variation  diminishing  property  of  G®  ensures  that  new 
contours  cannot  appear  at  positive  values  of  a,  nor  can  arches  be  closed  from  below. 

A  feature  is  delimited  by  right  and  left  inflection  points.  Each  closed  arch  of  a 
fingerprint  corresponds  to  a  feature.  For  example,  an  arch  composed  of  a  dotted  left 
contour  and  a  solid  right  contour  represents  a  hill.  Often  fingerprint  arches  are  nested:  a 
hillock  within  a  valley  would  give  rise  to  nested  arches,  the  hillock  corresponding  to  the 
inner,  lower  arch.  The  use  of  inflection  points  to  delimit  features  seems  more  natural  than 
using  relative  maxima  or  minima  of  elevation  to  delimit  hills  or  valleys.  At  each  sn^jothing 
scale  a,  the  inflection  points  of  F(x  ,a)  segment  the  x-axis  into  mutually  disjoint  hills  or 
valleys.  If,  instead,  we  were  to  consider  a  hill  as  delimited  by  adjacent  relative  minima,  and 
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a  valley  as  delimited  by  adjacent  relative  maxima,  then  hills  and  valleys  would  contain 
substantial  intervals  of  overlap. 

Sign  changes  of  the  inflection  points  in  the  upper  and  lower  parte  of  Fig.  1  have 
been  encoded  by  solid  and  dotted  lines.  A  closed  arch  formed  by  a  solid-dotted  line  is  a 
subtraction  feature  or  dip,  while  a  closed  arch  formed  by  a  dotted-solid  line  pair  is  an 
addition  feature,  or  tocal  peak.  Note  how  solid  and  dotted  lines  alternate  at  any  given 
value  of  sigma.  By  reading  the  fingerprint  from  higher  to  lower  values  of  sigma,  we  can 
describe  the  shape  and  properties  of  the  elevation  profile  with  an  increasing  level  of 
detail.  This  qualitative  description  of  the  elevation  profile  from  the  fingerprint  can  be  made 
much  more  quantitative;  in  fact,  it  is  possible  to  reconstruct  the  entire  elevation  profile 
from  the  fingerprint  to  within  a  scaling  parameter  and  a  linear  shift  (Yuille  and  Poggio,  Ref. 
3).  The  fingerprint  provides  a  method  to  identify  and  rank  key  features  of  an  elevation 
profile  and  a  formalism  to  describe  a  profile  to  arbitrary  levels  of  detail. 

III.  Summary  of  Program  Results 

General  properties  of  scale  space  analysis 

Initial  effort  on  this  program  focused  on  development  of  general  properties  of  scale 
space  analysis  that  are  independent  of  application  to  particular  forms  of  input  data: 

(a)  software  development  and  testing 

(b)  analysis  of  the  influence  on  the  fingerprint  of  software  parameters,  such  as  choice  of 
data  extension  method,  Gaussian  mask  size  or  mask  cropping 

(c)  formulation  and  testing  of  feature  properties  and  relationships  that  can  be  deduced 
from  the  fingerprint. 

Our  scale  space/fingerprint  software  was  designed  to  handle  very  general  input  data, 
so  that  the  analyses  (b)  and  (c)  were  based  on  application  of  the  software  to  a  variety  of 
input  data,  including  synthetic  curves  representing  common  edge  and  feature  types. 
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spectral  reflectance  curves  of  soils  and  common  terrain  materials,  terrain  elevation 
profiles,  and  spectroscopic  data  of  various  pathogenic  and  non-pathogenic  chemicals. 
General  scale  space  properties  developed  on  this  project  are  contained  in  publications  [1 , 
2,  3],  which  we  summarize  individually.  Abstracts  of  each  of  these  papers  are  included  in 
Section  V  of  this  report. 

1 .  M.  A.  Piech  and  K.  R.  Piech.  Symbolic  representation  of  hyperspectral  data.  Applied 
Optics  26  (1987),  4018-4026. 

The  fingerprint  provides  a  context  within  which  inflection  points  are  tracked  across  a 
continuous  range  of  filtering  scales.  Lesser  features  die  out  before  stronger  ones.  We 
can  spatially  locate  a  feature  by  the  left  and  right  locations  of  the  base  of  its  fingerprint 
arch.  The  scale  at  which  an  arch  closes  is  a  measure  of  the  magnitude  or  persistence  of 
the  feature.  This  measure  enables  us  to  compare  magnitudes  of  distant,  as  well  as  near, 
features.  Each  feature  can  be  described  by  a  triplet  consisting  of  the  scale  at  which  the 
teatuie  arui  closes,  and  ine  left  and  right  inflection  points  delimiting  the  feature.  This 
description  is  compact,  quantitative  and  hierarchical,  describing  the  hyperspectral  curve 
by  its  most  important  structural  features  first,  followed  by  features  of  lesser  importance. 

For  an  isolated  feature,  we  develop  a  scale/area  relationship;  the  scale  at  which  its 
corresponding  arch  closes  is  proportional  to  the  cube  root  of  the  area  enclosed  within  the 
feature.  Here  the  area  enclosed  within  the  feature  refers  to  the  area  between  the  graph  of 
f  and  the  straight  line  segment  connecting  the  inflection  points  that  delimit  the  feature. 
The  fingerprint  is  independent  of  measuring  units  chosen  for  the  input  function  f, 
although  the  choice  of  units  does  affect  the  constant  of  proportionality  of  the  scale/area 
relationship.  Moreover,  the  fingerprint  is  unchanged  if  a  linear  function  is  added  to  f,  so 
that  superposfiion  of  terrain  on  a  tilted  plane  leaves  the  terrain  fingerprint  unchanges. 

Three  methods  of  data  extension  are  compared;  (i)  end  point  value  extension,  (ii) 
average  value  extension,  and  (iii)  linear  fit  to  the  final  five  data  points  at  each  end  of  the 
elevation  profile.  We  argue  that  extension  by  linear  fit  is  more  appropriate  lor  fingerprint 


analysis,  since  only  this  method  does  not  introduce  spurious  inflection  points.  For 
numerical  convolution  with  a  Gaussian  mask,  the  mask  must  be  cropped  to  a  finite  width. 
The  mask  size  should  change  dynamically  with  the  standard  deviation  o  of  the  Gaussian 
mask,  and  should  be  no  less  than  3o  in  halfwidth.  Mask  cropping  is  shown  to  affect  the 
fingerprint  by  causing  a  drift  in  the  location  of  inflection  points  at  very  small  a.  especially 
where  the  input  curve  is  locally  linear.  As  a  consequence  we  choose  to  locate  the  left  and 
right  inflection  points  of  features  at  their  2a  positions. 


2.  M.  A.  Piech.  Comments  on  fingerprints  of  two  dimensional  edge  models.  Computer 
Vision,  Graphics  and  Image  Processing  A2  (1988),  381-386. 

Scale  space  images  and  fingerprints  of  the  most  common  two  dimensional  edge 
models  (step  edge,  pulse  edge  and  staircase  edges)  can  readily  be  deduced 
geometrically  from  fingerprints  of  the  corresponding  one  dimensional  edge  models. 

At  smoothing  scale  a,  the  range  of  influence  of  a  feature  edge  is  approximately  3o 
units.  A  practical  consequence  of  this  3a  effect  is  that  fingerprints  of  adjacent  edges  will 
be  identical  to  those  of  individual  isolated  edges  up  to  a  scale  of  one-third  the  distance 
between  the  edges. 

3.  M.  A.  Piech  and  K.  R.  Piech.  Hyperspectra!  interactions-  invariance  and  scaling. 
Applied  Optics  23  (1989),  481-489. 

Dominant  scale  space  features  (i.e.  pairings  of  inflection  points  that  persist  over  large 
smoothing  scale  ranges)  can  be  used  to  segment  a  profile  into  independent  regions.  The 
strength  and  location  of  features  in  one  region  (as  measured  by  fingerprint  analysis)  are 
independent  of  features  in  other  regions.  Consequently,  fingerprint  analysis  can  be 
carried  out  on  the  individual  profile  regions,  and  results  compared  across  the  regions. 

We  extend  our  previous  work  [1]  relating  scale  to  feature  area.  The  scale  at  which  a 
fingerprint  feature  vanishes  is  independent  of  feature  shape,  and  depends  only  on  the 
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area  of  the  feature.  Scaling  and  invariance  rules  are  developed  for  common  cases  of 
interacting  features.  At  large  feature  separations  the  fingerprint  beharior  of  one  feature  is 
independent  of  the  other  feature.  At  smaller  separations  two  features  interact  and  their 
scales  are  damped.  Below  a  critical  separation  distance  features  nest,  and  a  plot  of  feature 
scales  from  large  to  small  separations  exhibits  a  bifurcation  behavior.  Scales  of  features 
above  the  bifurcation  point,  scales  of  the  nested  features,  and  the  location  cf  the 
bifurcation  point  depend  only  on  feature  areas  and  not  on  individual  feature  shape  nor  on 
shape-associated  parameters. 

Terrain  modeling  results 

During  the  second  phase  of  this  program  we  correlated  fingerprint  and  fractal 
analysis  of  elevation  profiles  from  Digital  Elevation  Model  data  The  principal  results  are 
discussed  in  publication  4. 

4  M  A.  Piech  and  K.  R.  Piech  Fingerprints  and  fractal  terrain.  Mathematical  Geology  22 
(1990),  457-485. 

In  this  paper  we  describe  the  application  of  fingerprint  techniques  to  elevation  data 
for  three  major  landforms:  arid  shale,  humkj  shale,  and  granite  (Fig.  2).  Granite  data  is 
drawn  from  classic  examples  of  fault  block  mountains  in  the  Sierra  Nevada  range  of 
California,  and  from  the  Sawtooth  mountains  of  central  Idaho,  which  were  formed  by  deep 
and  extensive  erosion  of  the  Idaho  Batholith.  This  selection  enables  us  to  analyze  not 
only  gross  distinctions  such  as  those  between  granite  and  shale,  but  also  finer 
distinctions  between  arid  and  humid  shale,  and  between  granites  of  different  geophysical 
origins. 

The  simplest  characterization  of  topographic  features  from  the  scale  space 
fingerprint  is  the  feature  density  function  N(a),  which  represents  the  number  of  tetrain 
features  that  persist  beyond  a  smoothing  scale  a  (Fig.  3).  The  feature  density  function 
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Fig.  2  Typical  arid  shale,  humid  shale,  and  granite  fingerprints.  The  fingerprint  features 
(arches)  of  arid  shale  are  densely  packed  with  pattern  uniformity  at  small  values  of  a. 
Fingerprint  features  of  humid  shale  are  not  as  densely  packed  as  arid  shale  (85%  as  many 
features  as  arid  shale),  and  lack  pattern  uniformity  at  small  scales.  The  feature  density  of 
granite  fingerprints  is  roughly  half  that  of  either  arid  or  humid  shale,  and  feature  spacing  is 
uneven  at  most  scales. 
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Fig.  3  (Above)  Feature  density  functions,  extracted  by  fingerprint  analysis,  for  arid  shale, 
humid  shale,  and  granite.  The  granite  feature  density  function  is  drawn  for  o  >  2.5,  since  on 
this  domain  CA  granite  and  ID  granite  have  a  common  feature  density  function. 

(Below)  Fractal  scaling  of  the  feature  density  functions  for  arid  and  humid  shale  is  confirmed 
by  the  log-log  linearity  of  the  feature  density  functions  over  the  range  o  =  1  to  5. 


10 


differs  for  shale  and  granite  landforms  and  between  individual  expressions  of  the 
landforms  such  as  arid  and  humid  shales.  For  a<7  the  feature  density  functions  of 
granite,  arid  shale  and  humid  shale  are  clearly  ordered,  with  N(a)granite  ^  *^(‘^)humid  shale 
<  N(a)ang  shale-  f^urther,  the  [N(1 ),  N(1 .5)]  pair  of  the  feature  density  function  is  a  strong 
discriminant  betwen  shale  and  granite,  between  arid  and  humid  shale,  and  between  CA 
and  ID  granites.  The  ordering  of  the  feature  density  functions  and  the  number  of  features 
at  small  a  agree  with  intuitive  physical  concepts  of  the  prevalence  of  edges  in  the 
topography. 

There  is  an  intuitive  and  natural  correspondence  between  the  scale  changes  used 
to  establish  fractal  behavior  and  the  scale  space  smoothing  in  fingerprint  analysis.  Both 
concepts  evaluate  changes  in  behavior  as  a  function  of  measuring  or  smoothing  scale. 
We  therefore  attempted  to  evaluate  fractal  behavior  of  the  landforms  by  examining  the 
scale  space  fingerprints  and,  in  particular,  the  feature  density  functions  extracted  from  the 
fingerpints. 

We  first  applied  conventional  methods  of  fractal  analysis  to  determine  if  the 
landforms  could  be  described  by  fractal  parameters  over  a  range  of  physical  scales.  The 
conventional  methods  included  computation  of  trail  length  and  surface  area  as  a  function 
of  measuring  scale.  Both  the  length  and  area  computations  demonstrated  that  the  shales 
were  fractal  and  the  granites  were  not  fractal.  However,  these  conventional  methods  were 
surprisingly  insensitive  to  fractal  behavior  and  required  careful  treatment  of  the 
dimensions  of  coordinates  for  a  meaningful  description  of  fractal  behavior. 

The  fingerprint  modeling  results  in  a  new  expression  of  fractal  scaling  from  the 
feature  density  function.  There  is  strong  consistency  between  the  expression  of  scaling 
from  the  feature  density  function  and  expressions  of  scaling  by  fractal  parameters 
extracted  with  conventional  methods  such  as  area  scaling  and  trail  length  scaling.  The 
conventional  methods  and  the  feature  density  expression  agree  in  predicting  a  lack  of 
fractal  behavior  for  granite  landforms  and  in  predicting  fractal  behavior  for  shale  landforms. 
The  feature  density  expression  also  agrees  quantitatively  with  the  conventional  methods 
for  the  range  of  scaling  behavior  of  the  shale  landforms.  The  conventional  computations 
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of  scaling  behavior  result  in  scaling  for  measuring  yardsticks  up  to  180m  for  both  arid  and 
humid  shale.  The  feature  density  function  exhibits  scaling  behavior  to  a  o  of  5  units,  or 
150m.  Although  we  cannot  equate  a  measuring  scale  with  a  smoothing  scale,  because 
the  smoothing  scale  represents  a  standard  deviation,  these  scaling  ranges  are 
compatible.  The  fingerprint  approach  is  also  much  more  sensitive  to  fractal  scaling  than 
conventional  methods,  and  automatically  provides  consistency  in  scaling  units,  avoiding 
difficulties  encountered  in  conventional  methods. 

Extension  of  fingerprint  theory 

Scale  space  fingerp'int  theory,  as  employed  in  this  project,  operates  on 
one-dimensional  elevation  profiles.  Although  the  analysis  is  quite  stable  for  contiguous 
profiles,  it  would  be  very  useful  if  we  could  develop  a  fully  two-dimensional  (2-D) 
fingerprint  theory  which  would  enable  tracking  from  coarse  to  fine  scales  in  order  to 
accurately  locate  feature  boundaries.  Toward  the  end  of  this  project  we  studied  2-D  edge 
detection  theories  to  determine  whether  current  theories  could  serve  as  a  basis  for  a 
fingerprint  theory.  Publication  5  follows  from  these  studies. 

5.  M  .  A.  Piech.  Decomposing  the  Laplacian.  Pattern  Analysis  and  Machine  Intelligence, 
to  appear. 

In  a  recent  paper  on  authenticating  edges  detected  by  zero-crossings  of  second 

order  differential  operators,  Clark  [Ref.l]  derived  a  formula  which  shows  that  the  Laplacian 

V  f  differs  from  the  second  directional  derivative  ^  in  the  gradient  direction  of  f  by 

an2 

the  product  of  the  curvature  k  of  the  level  curve  for  f  through  (x,  y)  with  the  gradient 
magnitude  |Vf|.  In  this  note  we  present  an  elementary  proof  of  Clark's  formula. 

We  believe  that  an  entirely  new  approach  to  2-D  feature  detection,  based  on  the 
tracking  of  curvature  regions  with  increased  Gaussian  filtering,  can  serve  as  a  basis  for  a 
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truly  two-dimensional  scale  space  fingerprint  theory.  We  have  been  exploring  the 
behavior  of  zero  curves  ot  the  Gaussian  curvature,  as  the  surface  undergoes  progressive 
smoothing.  The  graph  of  these  zero  curves  versus  smoothing  scale  resembles  a  glove, 
and  we  would  hope  to  eventually  be  able  to  characterize  terrain  type  by  feature  density 
functions  (drawn  from  the  glove)  which  enumerate  the  areal  density  of  features  as  a 
function  of  their  persistence  across  smoothing  scales.  Such  a  2-D  fingerprint  theory 
would  be  a  more  appropriate  and  valuable  descriptor  for  terrains. 
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V.  Abstracts  of  Publications  Listed  Above 


Symbolic  Representation  of  Hyperspectral  Data 

We  devebp  a  symbolic  representation  of  hyperspectral  data  using  the  scale  space  techniques 
of  Witkin.  We  create  a  scale  space  image  of  hyperspectral  data  from  convolution  with  Gaussian 
masks  and  then  a  'fingerprint'  that  extracts  individual  features  from  the  original  data.  The 
fingerprint  provides  a  context  that  pairs  inflection  points  and  assigns  them  to  a  feature;  generates 
a  measure  of  importance  for  each  feature;  and  relates  features  to  each  other.  The  representation  is 
an  ordered  sequence  of  triplets  containing;  a  measure  of  importance  related  to  the  area  of  each 
feature,  and  the  left  and  right  inflection  points  of  the  feature.  The  description  is  compact, 
quantitative  and  hierarchical,  describing  the  hyperspectral  curve  by  its  most  important  structural 
features  first,  followed  by  features  of  lesser  importance. 


Comments  on  Fingerprints  of  Two  Dimensional  Edge  Models 

Shah,  Sood  and  Jain  have  recently  published  an  interesting  scale  space  analysis  of  pulse  and 
staircase  edge  models,  both  one  and  two  dimensional.  In  this  note  we  comment  upon  Shah,  Sood 
and  Jain's  analysis  of  the  two  dimensional  step  edge,  pulse  edge  and  staircase  edge  models.  Their 
derivation  of  the  scale  space  images  and  fingerprints  can  be  simplified  by  taking  advantage  of  a 
key  geometric  feature  of  Gaussian  filters,  namely,  rotational  invariance.  The  fingerprints  of  these 
three  models  can  easily  and  directly  be  geometrically  deduced  from  the  fingerprints  of  the  one 
dimensional  models.  The  fingerprints  should  be  viewed  as  cylinders  over  a  base  curve  which  is 
precisely  the  fingerprint  of  the  corresponding  one  dimensional  edge  model.  In  this  way  fingerprints 
of  the  two  dimensional  models  can  be  immediately  visualized  from  their  one  dimensional 
counterparts.  We  also  demonstrate  that  the  range  of  influence  of  one  edge  upon  another  edge 
located  a  distance  d  away  begins  at  a  scale  of  d/3. 


Hyperspectral  Interactions:  Invariance  and  Scaling 

We  present  new  invariance  and  scaling  results  for  scale  space  analysis  of  hyperspectral  data. 
First,  we  note  that  a  hyperspectral  curve  can  be  segmented  into  independent  regions  selected  by 
features  of  scale  space  fingerprints.  These  fingerprint  features  are  persistent  inflection  points  that 
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precisely  locate  major  atmospheric  features  that  define  the  regions.  The  strength  and  location  of 
hyperspectral  features  in  one  atmospheric  region  are  independent  of  features  in  other  regions;  as  a 
result,  hyperspectral  analysis  can  be  simplified  to  a  region-by-region  analysis.  We  than  generate 
simple  scaling  and  invariance  rules  for  features  within  such  a  spectral  region.  We  show  that  the 
scale  of  individual  features  is  independent  of  the  details  of  feature  shape  and  depends  only  on  the 
area  of  the  feature.  Interacting  features  in  turn  exhibit  a  fascinating  bifurcation  behavior:  at  large 
separations  features  behave  independently;  at  smaller  separations  features  interact  and  their 
scales  are  damped;  bebw  a  critical  separation  distance  (the  bifurcation  point)  the  features  nest. 
The  scales  of  features  above  the  bifurcation  point,  the  scales  of  the  nested  features,  and  the 
location  of  the  bifurcation  point  depend  only  on  the  feature  areas  and  not  on  shape-associated 
parameters  of  the  individual  features. 


Fingerprints  and  Fractal  Terrain 

Scale  space  fingerprints  provide  a  natural  formalism  for  modeling  terrains  because  they 
describe  a  hierarchy  of  edges  over  a  broad  range  of  physical  scales.  The  simplest  fingerprint 
descriptor  is  the  feature  density  function  N(<j) — ^the  number  of  terrain  edges  remaining  after 
smoothing  a  terrain  profile  with  a  gaussian  of  standard  deviation  a.  Feature  density  functions 
computed  from  Digital  Elevation  Model  data  are  not  only  clearly  different  for  landforms  such  as 
shales  and  granite,  but  also  for  landform  expressions  such  as  arid  and  humid  shale. 

Feature  density  also  provides  a  description  of  fractal  behavior  that  agrees  with  conventional 
fractal  computations  such  as  area  and  trail  length  scaling;  shales  exhibit  fractal  behavior  over 
scales  up  to  about  200  meters,  while  granites  are  not  fractal  over  any  meaningful  range.  The 
feature  density  evaluation  of  fractal  behavior,  however,  is  superior  to  the  conventbnal  methods:  it 
is  more  sensitive  to  the  expression  of  fractal  scaling  and  does  not  exhibit  a  dependence  on 
measuring  units  that  can  mask  fractal  behavior  in  the  conventbnal  methods. 


Decomposing  the  Laplaclan 

We  correct  and  give  an  elementary  derivation  of  Clark’s  decompositbn  of  the  Laplacian 
V^f(x,y)  into  a  second  directbnal  derivative  in  the  gradient  directbn  of  f,  and  a  product  of  gradient 
magnitude  by  curvature  of  the  level  curve  through  f(x,y). 
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